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This paper concerns approximation in the Chebyshev, or minimax 
sense such that (i) a minimax approximation implies a maximum number 
of zero error points separated by equal error extrema, and (ii) the approx- 
imating function can be so formulated that the disposable parameters 
are all the coefficients in a polynomial, which may however be part of a 
more complicated function the rest of ivhich is prescribed. Weighted minimax 
polynomial approximations can be included, by multiplying the approx- 
imated and approximating functions by the weight factor. Analytic methods 
are described which yield approximately equal error extrema. They are 
sufficiently simple so that they may sometimes compete with currently 
used iterative numerical methods, especially when the degree of the dis- 
posable polynomial is large. Their most probable utility concerns explora- 
tions of available accuracies over wide ranges of design parameters such 
as degree of disposable polynomial, interval of approximation, and coef- 
ficients in prescribed parts of the approximating function. 

I. INTRODUCTION 

This paper concerns approximation in the Chebyshev sense, over a 
prescribed interval x a 5S x ^ x b of a continuous real variable x. As 
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defined, an approximation in the Chebyshev sense is a minimax approx- 
imation — one in which the maximum error is as small as is possible 
within given constraints on the approximating function. Minimax 
approximations in which errors are weighted by a prescribed function 
of the independent variable can also be treated as Chebyshev approx- 
imations, by multiplying the approximated and approximating functions 
by the weight function. 

Frequently, but not always, approximation in the Chebyshev sense 
implies an error of the "equal ripple" sort illustrated in Fig. 1 — that is, 
a sequence of equal positive and negative extrema with monotonic 
variations in between. General necessary and sufficient conditions for 
this are not known. However, the following conditions are sufficient: 
the p disposable parameters of the approximating function are to be 
such that the approximation error can be made zero at p arbitrary 
points within the approximating interval. Referring to Fig. 1, the 
arbitrary error points divide the approximation interval into p + 1 
segments. There is to be a particular division such that the error function 
achieves its maximum magnitude p + 1 times — at the two edges of 
the approximation interval and once within each of the p — 1 interior 
segments. There are to be no other local extrema within the approx- 
imation interval. Generally, shrinking any one of the p + 1 segments 
(by bringing two zero error points closer together or one closer to an 
edge of the approximation interval) tends to reduce the corresponding 
error extremum. Conditions are to be such that all the p + 1 equal 
extrema can be reduced simultaneously only by shrinking all the 
p + 1 segments, which is impossible without shrinking the given approx- 
imation interval. These conditions are encountered in many practical 
problems and are assumed here. Thus we are concerned only with equal 
ripple approximations like Fig. 1. 

Exactly equal ripple approximations have long been known for a 
very few special cases (which have been useful for example in filter 
design). Iterative numerical methods have been developed for the 
solution of various more general problems and are described in text- 




Fig. 1 — An equal ripple error function (p zeros; p + 1 segments; p + 1 extrema). 



MI. MM AX APPROXIMATION 6 

books such as Ref. 1. In contrast, this paper describes analytical pro- 
cedures which yield error extrema of approximately equal amplitude. 
Their full range of validity has not been determined. However, they 
are clearly appropriate for a substantial, although poorly defined class 
of problems. It is characterized further later. 

Useful applications are likely to concern equal ripple problems which 
have not been solved exactly by analysis and which involve so many 
disposable parameters that iterative numerical solutions are likely to 
be more costly. The most useful applications probably concern prelim- 
inary explorations over primary design parameters (such as intervals 
of approximation, magnitudes of errors, and degrees of approximating 
functions) before numerical refinement of specific designs. Accordingly, 
this paper emphasizes relatively simple means for approximating equal 
ripples and says little about more complicated higher order approx- 
imations. 

The procedures apply only to approximating functions characterized 
as follows. The disposable parameters must be all the coefficients in a 
polynomial (which may have been obtained, however, by some sort of 
transformation on the original independent variable and/or the approx- 
imating function). This is referred to as the disposable polynomial. 
On the other hand, the disposable polynomial may be only a part of a 
more general approximating function the rest of which is prescribed 
in advance (for example, the numerator of a rational fraction with a 
prescribed denominator). Weighted as well as unweighted minimax 
approximations are included. For some problems, closed form formulas 
are obtained for approximate error size as functions of the degree of 
the disposable polynomial, usable for degrees of any size. For other 
problems, the error size is related to an eigenvalue of a certain matrix 
equation, but the order of the matrix may be small even though the 
degree of the disposable polynomial is arbitrarily large. 

A primary concern here is the distinction between simple truncation 
of infinite series of Chebyshev polynomials and approximation in the 
Chebyshev or minimax sense. The functions which we are to approx- 
imate can be expanded into infinite series of Chebyshev polynomials. 
Approximations with polynomials of degree n can be obtained by 
simply truncating the infinite series after the terms of degree n. How- 
ever, simple truncation does not usually give an approximation of the 
minimax sort. A polynomial of degree n which approximates the given 
function in the minimax manner can be represented as a linear combi- 
nation of Chebyshev polynomials, but the coefficients are usually 
different from those in the truncated infinite series. 

One way to approach approximation in the Chebyshev sense is to 
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start with the truncated series of Chebyshev polynomials. Then cor- 
rections to the coefficients are determined, to obtain equal ripple error 
functions. Such a procedure has been used before, for example in Refs. 
2 through 4, and is used here. Departures from the previous work 
known to the author include simple approximations to ideal solutions 
formulated for more general approximating functions and for weighted 
as well as unweighted minimax approximations, as opposed to more 
rigorous analyses of more restricted problems. 

Sometimes truncation of an infinite series of Chebyshev polynomials 
yields an approximately equal ripple error function without further 
adjustment of the coefficients. The procedures for adjusting the coef- 
ficients, described herein, sometimes also give an initial insight into 
whether or not adjustments are needed. 

It is interesting to note that some 35 years ago a conference was 
held in the office of T. C. Fry, at Bell Telephone Laboratories, to 
consider some filter patents offered for sale by W. Cauer. One of the 
patents disclosed Cauer's equal ripple image impedance and transfer 
functions, which soon became famous among circuit theorists, but did 
not include proofs or derivations. At the conference, S. A. Schelkunoff 
asserted a very simple principle which enabled him to confirm and 
interpret Cauer's formulas. However, it did not explain how Cauer 
might have derived or discovered the formulas. The principle applies 
also to more general equal ripple approximations. It does not, by 
itself, solve the approximation problem, but it does furnish a starting 
point from which to develop procedures which do. We call it Schelkunoff's 
principle. 

Section II describes Schelkunoff's principle. Section III solves two 
problems for which exactly equal ripple solutions are easily found. 
Section IV develops general procedures, whereby approximate so- 
lutions can be obtained for a large class of problems. Section V further 
clarifies the general procedures by means of examples. 

Various aspects of the procedures described here bear some relation 
to other work. Section VI notes some of these relationships. Finally, 
Section VII reviews and summarizes the general conclusions, including 
a comment on the possibility of generalizations to disposable rational 
fractions. 

ii. schelkunoff's principle 

Consider first a function T n (x) proportional to a Chebyshev poly- 
nomial, defined by 
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T„(.r) = r^ri cos (n cos ' x) . 



(1) 



It is illustrated in Fig. 2a, for n = 4. It may be regarded as an "equal 
ripple" approximation to zero, over the interval —1 ^ x fS +1, by 
a polynomial of degree n in which the coefficient of x n is required to 
be K„ . Let 



X = COS cp. 

Substitution in equation (1) gives 

7' 

T n (.r) = T n (fp) = 2^i cos n<p. 



(2) 



(3) 



The new function is illustrated in Fig. 2b, again for n = 4. Note that x 
is periodic in <p with period 2t and T n (<p) is periodic in <p with period 
2w/n. Thus there are n periods of T n (ip) in each period of x. 

Stated with a little more detail, we have this situation: The original 
function T„(:r) has "equal ripples" in the sense of equal extrema. How- 
ever, the extrema are not uniformly spaced and hence the ripples 
differ as to width. The periodic transformation from x to <p has two 
important properties. As <p increases, x sweeps back and forth across 
the approximation interval, —1 ^ x ^ +1. In each interval in which 
x varies monotonically from ± 1 to =F 1 the <p scale is a distortion of the 
x scale such that the ripples of T n (ip) are uniformly spaced and are 




Fig. 2 — Illustrating Schelkunoff's principle. 
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also otherwise identical. The ripples could not be made identical bj- 
any distortion of the x scale alone if T n (x) had unequal maxima or 
unequal minima. This is a special case of Schelkunoff's principle. 

More generally, let E(x) be a function of x with the following prop- 
erties over an interval x a ^ x ^ x b : The function E(x) is real and single 
valued; there are a number of local maxima, all equal; there are a 
number of local minima, all equal; the equal extrema for the interval 
include the end points E(x a ), E(x b ) (at which dE/dx need not =0).* 
Then Schelkunoff's principle asserts the existence of a transformation 

x = r(«0 (4) 

with the following properties: The original variable x is periodic in the 
new variable <p; as tp increases x sweeps back and forth over the given 
interval x a ^ x ^ x b , monotonically each way once each period; the 
periodic function 

E(<p) = E[r(*)] (5) 

has a number of periods in each period of x, equal to one less than the 
number of extrema of E(x) in the given interval of x (including the end 
points). In applications to approximation in the Chebyshev sense, 
E(x) and E(<p) represent the equal ripple error, as functions of x and <p. 

The transformation x = T(ip) clearly is not unique, for there are 
obvious transformations on <p itself which retain the desired character 
of E(<p). For example, tp can be replaced by <p + q((f), where q(<p) is 
periodic with the same period as E{tp) and is such that <p + q(<p) is 
monotonic in tp. When E(x) is continous (in the given interval of x), 
a particular <p + q(<p) will make E(<p) sinusoidal. 

We do not attempt a very general, rigorous proof of Schelkunoff's 
principle; we merely use it as a guide to a strategy for solving minimax 
problems. However, a demonstration of the principle for a specific 
class of problems will be implicit in what follows, for we shall find 
transformations which do in fact change our equal ripple errors into 
sinusoidal errors. 

In later sections we will again use the transformation (2), or a gen- 
eralization for end points other than x = ±1. Usually, however, it 
will not be a Schelkunoff transformation. We will use it to transform 
the disposable polynomial in x into a finite Fourier series in <p. The coef- 



* Problems can bo found such that minimax approximations have equal ex- 
trema which do not include both end points. Then the number of local extrema 
for the approximation interval is abnormally large when the end points are 
counted. Such problems arc not considered here. 
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ficients of the Fourier series are to be chosen in such a way that the 
overall approximation is approximately sinusoidal on a distortion of 
the tp scale. Similar strategies have been used before, for example in 
Refs. 2 through 4. 

Means for determining the distortion of the <p scale and the adjustment 
of the Fourier coefficients are introduced by means of two examples in 
the next section. 

III. TWO GENERALIZATIONS OF CHEBYSHEV POLYNOMIALS 

The two problems described below are solved exactly. The form of 
the solutions suggests approximate solutions to more general problems. 

3.1 A Rational Function Generalization of Chebyshev Polynomials 

Consider the following generalization of Chebyshev polynomials: Let 

TM - g; (6) 

in which P(.r) is a polynomial of degree n and T>(x) is a polynomial of 
degree ^ n. Suppose D(x) is prescribed in advance and P(.r) is to be 
chosen in such a way that T 0n (s) has equal ripples like those of a 
Chebyshev polynomial in the interval -1 ^ x ^ +1. More specif- 
ically, require that T Dn (x) = ±1 at n - 1 local extrema within the 
interval -1 ^ x ^ +1 and at the end points x = ±1. Real zeros of 
D(x) are to be excluded from the interval -1 ^ x ^ +1. The con- 
ditions on the extrema insure that all n of the zeros of P(.r-) will be in 
the interval. 

Let tp be defined again by equation (2) and note that the real axis 
in the tp plane corresponds to the real interval —1 ^ x ^ +1 in the 
x plane. If cos tp is written in terms of exponentials, the polynomial 
P(.t) can be related to <p by 

P(.t) = P(e i,p ) + P(c~'*) (7) 

in which P(-) is a polynomial of the same degree, n, as P(-)- Given the 
coefficients of P(-) it is a simple matter to compute the coefficients of 
P(-). We shall consider our problem solved when we have found the 
coefficients of P(-) required for our equal ripple conditions. 

It is convenient to relate the prescribed denominator D(.r) to <p in 
the following slightly different way: 

D(.r) = D(e") D(f iv ) 
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in which D(-) is a polynomial of the same degree, ^ n, as D(-). If 
D(-) and D(-) are written in factored form there is a one to one corre- 
spondence between factors. Thus if (x„ — x) is a factor of D(z) and 
(1 — y a e" ( ') a corresponding factor of D(e ,lp ), 

x a - x = M a {\ - 7 „0(1 - y* tr ) (9) 

in which M a is a constant scale factor. 

By equation (2), e 1 * 1 = ±1 at x = ±1, and hence 

x a ± 1 = M a (l ± y,) 2 . (10) 

Given x„ , two solutions for y a are easily obtained, for which | y„ | is 
respectively < 1 and > 1. (Exclusion of zeros x a of D(x) from the real 
interval — 1 ^ x ^ +1 removes the possibility of | y a | = 1.) We need 
the solution for which [ y„ | < 1, for reasons which will soon be apparent. 
From equation (11), with the sign of the square root chosen for [ y„ | < 1, 



x a — 



\l -fertf>o- 



(11) 



1 + y„ x. + 

The scale factor M a need not concern us at this time. 

The function T Dn (x) can now be mapped into a function T Dn (<p) in 
terms of equations (7) and (8). 

T Dn{ x) = T M = P(ein+P(e ~ i *\ (12) 

D(e 1v ) Die"") 

By Schelkunoff's principle, our requirements on the extrema of T Dn (x) 
imply that T Dn (<p) has the following special form 

m P(,")+P(« : ") 

Die'*) D(e-'«) 
= i e <[»P+/( P )i _j_ y-i [»»+«»)]_ ( 13 ) 

The variable <p -f- (l/n)/(v) is a distortion of the <p scale for Schelkunoff's 
principle, for which f(<p) is to be periodic in <p with period 2t and <p + 
(l/n)f(cp) is to vary monotonically with <p. 

Given j{ip) one can easily find P(e i,p ) by equation (13). The problem 
is to find an f(<p) for which P(-) is a polynomial of degree n. From 
equation (13) 

P(e") + P(e-*) = J 2 * v n {e)v n {e ) I (w) 



If P(e ,> ) is to have no terms in e i " p with o- > n, e i/(r) needs to cancel 
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out D n (e iv ). This suggests 



e 



/<¥>> _ 



Z)(e-') 



(15) 



W) 

(and note that this does make f(<p) real when <p is real). Substitution 
in equation (14) gives 

P(e") + P(r'*) = §e' n * DV^) + K"** W). (16) 

Expanding the right side gives a polynomial in e M *. When polynomial 
Z)(*) is of degree ^ n (as assumed) there "wall be no powers of e"' out- 
side the range — n to -\-n. Collecting positive powers (and half the 
constant term) gives P(e'*). 

The function /(p) determined by equation (15) is periodic in ip with 
period 2w provided the zeros of D(X) he outside the unit circle in the 
X plane, which is assumed by equation (11). It is easily shown that the 
same condition makes <p + (l/n)/(<p) monotone in <p. Once P(- ) is known 
it is a simple matter to find P(x) by means of equations (2) and (7). 
It is probably simplest to omit the scale factor M, of equation (10) in 
the initial formulation. This does not affect the ratio in equation (15), 
but only the scale factor of the polynomial P(x), which can be corrected 
later on [for example to meet the condition T Dn (l) = 1]. 

Obvious generalizations of the problem include the following: For 
extrema A =b J (instead ± 1) use 

F(*) = A ± JT Dn (x)- (17) 

In a more general interval of x, say x a ^ x ^ x b , replace equation (2) by 

Xh "T X a . Xh X a /io\ 

x = ^ 1 g — cos * <> 18 ) 

and change equation (10) to 

x t - x a = M„(l - 7,) 2 , x a - x b = M.{\ + y,) 2 . (19) 

The function T?(x) defined by equation (17) has long been used by 
filter theorists, but previous derivations have been quite different. 5 
The form of equation (16) suggests a similar solution to the problem 
described below. 

3.2 An Irrational Generalization of Chebyshev Polynomials 
Now let 

t -w = S 1 (20) 



10 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1970 

in which P(z) is again a polynomial of degree n but S(x) is a poly- 
nomial of degree ^ 2n. Suppose S(x) is prescribed and that T Sn (x) is 
to meet the same conditions as to extrema as T Dn (x) in the previous 
subsection. In place of equation (8), we can now use 

S(.t) = S&^Sie-'*) 
to determine a polynomial £(•)• We can then replace equation (15) by 

e ,„ = \s(rnv (21) 

and then equation (16) by 

Pie") + Pie-'*) = y nv S(e- { n + K' n *<S(e"). (22) 

This makes P(-) again a polynomial of degree n. Note that T Sn (x) 
cannot be used in place of T Dn in equation (17), with A ^ 0, without 
changing the polynomial character of the numerator. 

IV. GENERAL FORMULATION'S 

This section shows how a large class of minimax approximations can 
be approximated by generalizing the manipulations described above. 
In Section V, we clarify the general procedures further by providing 
examples. 

4.1 Unweighted Minimax Approximations 
Let 

P(z) = F(x) + e(as), x a ^ x S x b (23) 

in which P(x-) is a disposable polynomial of degree n, F(x) is a given 
function to be approximated by P(x) in the interval x a ^ x ^ x b , and 
t(x) is the error in the approximation. For what P(x) is t(x) smallest in 
the minimax sense? We assume that the minimax t(x) has the equal 
ripple form (Fig. 1) and we seek only approximations to equal ripples. 
We also restrict the class of applicable functions by certain further 
assumptions which can best be introduced a little later. 

As before, let x and <p be related by equation (18), so that x„ £ X £ x b 
maps into real <p, and replace P(z) by 

P(.r) =P(6")+P(0. (24) 

If P(-) is again a poljmomial of degree n, it is uniquely determined by 
P(-) [and x a , x b in equation (18)]. Now, however, we find it expedient 



MINIM AX APPROXIMATION 11 

to permit P(z) to include negative powers of z, up to z~ n . Thus, in 
equation (24), 

P(0 = J2 PJ". (25) 

This P( • ) is not uniquely determined by P( • ) ■ However, P( • ) is uniquely 
determined by P(-), and it is still a polynomial of degree n. We solve 
the approximation problem by finding a suitable P(-)> from which 
P(-) can be easily determined. 

We require that the mapping from x to p maps F(x) into a function 
of <p with a convergent Fourier series. This amounts to requiring that 
F(x) can be expanded into a convergent series of Chebyshev poly- 
nomials (denned to fit the given interval of approximation). Because 
equation (18) is even in <p, the Fourier series has only cosine terms. 
Then, replacing cosines by sums of exponentials, 

F(x) = F(e iv ) + F(f ,r ) 

(26) 

W) = £ W 

ff = 

in which the series expansion of F(e ,lfi ) converges when <p is real. 
The desired equal ripple error can be written 

t(x) = e cos [(n + 1) 9 + /(*)] (27) 

in which /(<£>) is again periodic in <p and represents the distortion of the 
cp scale per Schelkunoff's principle. In an equivalent exponential form 

t(x) = E(e ip ) + E(e-**) 
E( e <v ) = | 6 "<-+"* + «*>\ (28) 

The exponent i(n + 1)^, instead of imp as in the previous section, 
reflects the following circumstances: If the Chebyshev polynomial 
series corresponding to equation (26) is truncated after the polynomial 
of degree n, the first omitted polynomial is of degree n + 1. If all the 
other omitted polynomials have sufficiently small coefficients, the 
truncation error will approximate E(.t) of equation (23) with j(ip) = 0. 
Note also that a disposable polynomial of degree n has n + 1 disposable 
coefficients. These are an example of the p disposable parameters in 
the more general description of equal ripple errors in Section I. 

Using equations (24), (26) and (2S) in equation (23) gives 

P(e") 4- P(a~") = F(e iv ) + F(f i9 ) + E{e i9 ) + E(e~ ir ). (29) 



12 THE BELL SYSTEM TECHNICAL J0TJENAL, JANUARY 1970 

We arbitrarily equate the terms in exp (+*V) separately, so that 

P(e iv ) = F(e ip ) +E(e ir ). (30) 

If equation (30) is satisfied at all real <p so is the corresponding equation 
in exp (—i<p). Thus a solution of equation (30) is a solution of equation 
(29). But the converse is not necessarily true. Frequently, an exactly 
equal ripple approximation corresponds to a solution of equation (29) 
which is not a solution of equation (30). However, we will find that 
approximations with approximately equal ripples can frequently be 
derived from equation (30) , and in a much simpler way. 

In equation (30), expand P(-), F(-), E(-) per equations (24), (26) 
and (28). The result can be rearranged as follows: 

X-l ^ X-0 

G x = P n+1 -x - C n+1 _x , X^n+lX^rc + 1; 

= P n+1 _ x , n+ 1 < X ^ 2n+ 1. (31) 

In this equation, C n +i+\ is fixed by equation (26) but P n+ i_x is a dis- 
posable parameter in equation (25). Thus we seek an e and exp [ifiv)] 
with the following properties: First, (e/2) exp [if(<p)] is to be expandable 
in terms of positive and negative powers of exp (i<p). Second, the coef- 
ficients of positive powers are to cancel the corresponding coefficients 
C n+ i +x in equation (31). Third, the coefficients of negative powers are 
to be such that, with an appropriate e, |exp [if (<p)]| = 1 when <p is real, 
so that /(<p) is real and the error extrema are equal per equation (27). 
Sometimes it turns out that there are no negative powers beyond 
— (2n + 1). Then the left side of equation (31) can be adjusted to match 
the right side. In many other problems, approximately equal error 
extrema can be obtained by simply ignoring terms in negative powers 
beyond — (2n +1). 

Now consider the class of functions P(-) such that, in equation (31) 

t C. M e a ' = fg3 (32) 

in which A(-) is a polynomial of degree m and B(-) is a polynomial 
of degree n. If the series converges, as assumed, the zeros of A(z) will 
he outside the unit circle. 

Under conditions which we shall examine further, the appropriate 
exp [if(<p)] is now as follows: 
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, /( „ = A(e-'* )X(e'*) • (33) 

Aie'^Xie-*) 

in which 5 is small (at real <p) and X(-) is a polynomial determined by 
two further conditions. First, the zeros of X(z) must lie outside the unit 
circle. Second e and X(-) must be such that 

2 A{e i9 )X(f iv ) A(e«) X(e _,> ) 

in which N(-) is a polynomial. Let us examine the implications first 
and the existence of such an e and X(-) thereafter. 

When ip is real exp (i<p) and exp (—i<p) are conjugates, and so are 
identical polynomials in these two variables. This makes /(<p) real in 
equation (33), except for small corrections due to 5. When the zeros 
of A(z) and X(z) lie outside the unit circle, as required, the unit circle in 
the z plane maps into contours in the polynomial planes which do not 
enclose 0. This makes /(^) periodic in <p. 

The condition on the zeros of X(z) also permits the right side of 
equation (34) to be expanded: 

e- ip N(fi- tp ) m ^ , (35) 

X(e-'*) 7=1 

Using equations (32), (33), (34) and (35) in equation (31) now gives 

«(e"'*) + E $**" = 2 E G*" r J (36) 

(f=i ir-i 

P a = & + C n+1 _, , a in+lj 
= G a , n + 1 < <r < 2n + 1 ; 

2n + 2 

This 5 is small provided the zeros z,- of X{z) are such that zJ i2n+2) is 
small. When 5 is small, the actual error extrema will differ from e, but 
by no more than ± |5| e. 
Equation (34) requires 

| A(f i9 )X{e i9 ) + B(e iv )X(e- iv ) = A(e ip )e- i9 N(e" i9 ) . (37) 

The appropriate degree 77 of polynomial X(-) turns out to be one less 
than the number of poles of zB(z)/A{z) (including any poles at z = w). 
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When the degree of A(-) is greater than the degree of B(-) and the 
zeros of A ( • ) are distinct, a set of 77 + 1 homogeneous equations in the 
coefficients <?,• of X(-) can be derived by evaluating equation (37) at 
the zeros of A ( • ) . Then 

(m a + I Al))Q = (38) 

in which Q is the column matrix of the coefficients g, of X and M A 
and M B are square matrices of order 77 + 1. Under other conditions 
an equation of the same form can be obtained in other ways. 

Equation (38) requires e/2 to be one of 77 + 1 eigenvalues for which 
the matrix coefficient of Q is singular. Each eigenvalue determines a 
polynomial X(-) [including an arbitrary scale factor which cancels 
out in equation (33)]. For our purposes, we must choose an e which 
is real and such that the zeros of X(z) lie outside the unit circle. This 
raises a question of the existence of a suitable e and X(-). 

When degree 77 = 0, X(-) is a constant, equation (38) is a real linear 
equation in e, and there is no zero of X(-). When 77 = 1, X(') is linear, 
e is a root of a quadratic equation. It is not hard to show that the two 
roots are real [under our assumptions regarding zeros of A(-)] and 
that one (the larger) yields a zero z x of X{z) such that \z^ > 1. Equality 
occurs only in singular cases such that n + 2 zero error points are 
possible [even though there are only n + 1 disposable coefficients of 
P(x)] and can be so placed that there are n + 3 equal error extrema, 
instead of only n + 2. This may be seen by assuming that the zero of 
a linear X(z) is ±1, and then noting that, in equation (33), 

MD. = i±C = ±e .>. (39 ) 

X(e~'n l±<f" v " 

Conditions for the existence of a suitable e have not been established 
for 77 > 2. They are probably at least closely related to the (unknown) 
general conditions under which the minimax approximation has the 
equal ripple form of Fig. 1. 

The procedures described here are appropriate only when a suitable 
e does in fact exist. However, degrees 77 = and 1, for which existence 
has been established, are sufficient for many practical problems. For 
approximately equal error extrema, equation (32) itself need be only 
an approximation, and polynomials .A(-) and B(-) for which 77 = 
or 1 are likely to give a good enough approximation. This is particularly 
true when degree n of P(.r) is sufficiently large so that coefficients 
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C a , <x > n, in equation (26) approach a simple asymptotic behavior. 
Percentage variations between error extrema need not have to be 
very small even though the absolute errors must be very small. For 
example, a 10 percent variation between very small extrema may be 
acceptable, compared with large variations obtained by truncation 
of the infinite Chebyshcv polynomial series. 

Table I indicates the degree m of A(-) and n of B(-) for which tj = 
or 1. The column headed in + m + 1 indicates the number of dispos- 
able parameters in the rational fraction A/B which can be adjusted 
to approximate the sum in equation (32). 

The procedures for tj = and 1 are particularly well suited for rapid 
explorations of available error magnitudes as functions of initial design 
parameters, such as degree of the disposable polynomial, extent of the 
approximation interval, and parameters in the approximated function. 
When only the error magnitude is needed, it is not necessary to calculate 
the coefficients of polynomial P(-), which requires the series expansion 
(35). When tj = 0, the error magnitude is (approximately) the single e 
determined by equation (38). Then simple closed form formulas can 
frequently be obtained (and will be included in 4 of the 5 examples 
in Section V). When tj = 1, e is one of the two roots of the quadratic 
equation required by equation (38). (To meet the condition on the 
zero of X(-), the larger e must be chosen.) 

When e has been determined, it can be compared with the error e r 
obtained by simply truncating the Chebyshev polynomial expansion 
of F(.r). In terms of equations (26) and (32) 



e T = 



B(c j n 



+ 



B(C " P ) -,(„+!) «■ 



(40) 



A(e iv ) " A(e-*) 

Comparing the maximum e T (at real <p) with e indicates the improve- 



Table I — Value of in and n for which tj is or 1. 



Degree m 


Degree n 


Degree'?; 




olA(-) 


ofB(-) 


of A'(-) 


TO + M + 1 











1 


1 








2 


2 





1 


3 





1 


1 


2 


1 


1 


1 


3 


2 


1 


1 


4 
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ment to be obtained by the minimax refinement of the truncated 
series. Frequently, max e T occurs at <p = or ir, and then 

maxer = 2 f(Si)- ( 4i > 

4.2 Weighted Minimax Approximations 
Let 

P(z) = F(.r) + ^ e(x) (42) 

in which P(#) is again a disposable polynomial of degree n, T?(x) is 
again a given function to be approximated in the interval x a ^ x ^ x b , 
and the new function W(x) is a given weight factor. For what P(z) 
is e(x) smallest in the minimax sense? We will again assume that the 
minimax e(x) has the equal ripple form and will seek only approxi- 
mations to equal ripples. We will also assume that W(x) is bounded 
and positive definite in the approximation interval. (A point where 
W(x) = or oo would probably spoil the equal ripple character of the 
minimax approximation.) 

Map from x to <p as before and define P(-), F(-) and E(-) again by 
equations (24), (26), and (28). Express W(x) also in terms of expo- 
nentials, but as a product instead of a sum. More specifically, let 

H(x) = -log W(x) = H(e i,p ) + #(e~ ,p ) (43) 

and assume that W(x) is sufficiently smooth, as well as bounded and 
positive definite, so that H(z) is regular at \z\ ^ 1. Then 

^ = WW-), (44) 

with log D(z) regular when \z\ ^ 1. This D(«) is a generalization of the 
D(-) of Subsection 3.1 and of [£(•)]* of Subsection 3.2. Frequently 
it can be found by direct factorization of a function of e'* as in Section 
III. 

Equations like (29) and (30) can now be obtained as before. The 
only difference is that E(-) must now be multiplied by the product 
of functions of ^ in equation (44). Then equation (30) becomes 

P(e ,> ) = F(e iv ) + D(e^)D(e- i,p )E(e iv ) (45) 

and equation (31) becomes 
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2n+l 



S G*"' = 5 D(e")i)(e-' V /<r) + S C-i+xc"* ; (46) 

X-l & X-0 

(?x = P„ +1 _x - C n+1 _x, X £n + l; 

= F n+1 _ x , n + 1 < X ^ 2n+ 1. 

Retain the rational fraction B(-)M(0 of equation (32), but change 
equation (33) to 

</{r> = Bfe-'W'W) ^ (4?) 

2)(c'*)il(6**)Z(e-'*) 

so that 

Dti^DUr'W' = ^V'^OT'TW*) + a(e * f) (48) 

A^'We"") 
in which 5 and 5 are small. Then change equation (34) to 

e flyMfr—W) B^l) = e-'^Nje-") (4Q) 

2 A(e , '*)Z(e _,> ) A(e l> ) Z(e _,> ) 

Using equations (32), (47) and (35) in equation (46) now gives equation 
(36) again. From equation (49), equation (37) must be changed to 

| If(fi- tr )A(e- <9 )X<fi ir ) + B{e i9 )X{e- iv ) = A^e'^Nif 1 *). (50) 

Equation (50) can be used to find e, and the X(-) and N(-) needed 
for equations (35) and (36) . 

4.3 More General Approximating Functions 
Let 

*\P(x), x] = G(x) + t(x) (51) 

in which \?[P(.r), x] is a given function of x and a disposable polynomial 
P(x), G(x) is a given function to be approximated by ^[P(x), x] in 
the interval x a ^ x ^ x b , and e(.t) is the error in the approximation. 
For what P(z) is t(x) smallest in the minimax sense? Under certain 
further assumptions regarding \1>(- , •) this approximation can be 
transformed into a weighted minimax polynomial approximation. 
Assume an inverse SI'" 1 of ^[P(x), x], with respect to P(z)> exists 
over the approximation interval. Then equation (51) can be replaced 
by 

P(«) = ST 1 {[G(s) +•(*)],*}. (52) 
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Assume ^ f_1 (- , •) is sufficiently smooth and t(x) sufficiently small to 
justify the following approximation (in the interval x a ^ x ^ x b ): 

P(x) = V-^Gix), x] + 89l ^' X] «(*)• (53) 

This is in the general form (42) with 

F(a:) = * _1 [G(:r),;r] 
1 dV^lGix), x] 



(54) 



W(x) dG(.r) 

Thus Subsection 4.2 can now be applied provided the F(-) and W(-) 
determined by equation (54) meet the appropriate conditions. Recall 
that we required W(x) to be bounded and positive definite over the 
interval of approximation. However, reversing the sign of W(x) merely 
reverses the sign of t(x). Hence, in equation (54), we need only require 
that the partial derivative must be bounded and either positive definite 
or negative definite and sufficiently smooth for log D{z) to be regular 
when \z\ < 1. 

As a first example of the inversion of equation (51), let 

W(z)P(.t) = G(x) + t(x) (55) 

where W(x) and G(x) are given functions of x. Then 

= G(x) + e(.r) = G(x) _1_ 
* {X) W(x) W(.r) + W(») {h { ; 

As a second example, let 

[A(x) + B(z)P(z)]* = G(aO + t(x) (57) 

where A (re), B(rc), and G(x) are given functions of x, with B(x) and 
G(x) positive definite over the approximation interval. Solving for 
P(x) gives 

w-**te±a + .g*9 + ^». (58) 

If the term in e 2 (x) is omitted 

[A(x) + B(*)P(*)]« ^ G(.r) + e(.r) - ^ (59) 

in which terms in t"(x) have been neglected for a > 2. An equal ripple 
c(a;) in equation (57) yields an approximately equal ripple error if the 
last term is somewhat smaller than the extrema of e(x). 
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4.4 Relation to Phase Modulation 

The function of ip defined by equation (28) is similar to functions 
of time t used in communication theory to describe phase modulated 
signals. If <p is replaced by t in equation (28), 

E(e u ) « fl «»-«"+"<"i. (60) 

This is the exponential representation of a phase modulated signal 
in which the carrier (radian) frequency is n + 1 and the baseband 
signal f(t) is periodic with one period every n + 1 periods of the carrier. 
The signal may be regarded as the carrier plus sequences of upper and 
lower sidebands. The upper sidebands are determined by the coeffi- 
cients C n+1+X in equation (31). The lower sidebands are determined 
by the requirement of a purely phase modulated signal. Finally, if 
the sequence of lower sidebands extends as far as the negative carrier 
frequency — (n + 1) we truncate it at —n. 

Weighted minimax approximations can be interpreted similarly, 
in terms of simultaneous phase and amplitude modulation. 

4.5 Alternative Procedures 

It is obvious that the procedures described above can be varied in 
many different ways. A very few of the possible variations are noted 
below. 

Preliminary manipulations may be needed to obtain a formulation 
in which the disposable part is a polynomial. Also, the pertinent Fourier 
series may be sums of sines instead of cosines. Both these situations 
will be illustrated by Example 5, in Section V. 

If d^~ l [G(x), x]/dG(x) is expressed as a product of functions of x, 
D(e i,p ) can be formulated as a product of corresponding factors. A 
factor of the form (1 — x/x t ) p contributes a factor of the form 
[M a (l — y a e iip )] p , as in Section III. More generally, there may be 
advantages to replacing the D(-) of equation (44) by D(-) defined 

W ±2 (.r) = f>(e iv )f>(e- i,p ). (61) 

Then the D 2 (-) in equations (48), (49) and (50) is replaced by D ±l (-). 
It may sometimes be convenient to express P(x) and F(.r) as products 
of factors in exp (±*V) instead of sums, say 

P(x) = P(e i9 )P(fT"), F(x) = F{e iv )F(f iv ) (62) 

in which P(-) is a polynomial of degree n with no negative powered 
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terms. If a term in e is neglected, one can now replace equation (29) 

by 



P(e iv )P(e- iv ) = W) + ^P- U"") + ®te-^l . 
L F(e-")JL W)J 

Equating factors separately replaces equation (30) by 



(63) 



P(e <v ) = F(e iv ) + i^fe-?.. (64) 

V F(<f ,p ) 

Subsequent modifications of our previous procedures are now easily 
worked out. 

It would be possible to replace equation (33) by other functional 
forms for exp [if(<p)]. The moduli must approximate unity at real tp 
and expansions in positive and negative powers of exp (itp) must exist. 
Disposable parameters are to be adjusted so as to approximate the 
required coefficients of positive powers. However, except for very 
special functional forms [such as equation (33)], the adjustment is 
likely to be a quite complicated task. 

V. EXAMPLES 

This section further clarifies the general procedures by means of 
five examples. 

5.1 Example 1 
Let 

P(») 



(1 - */*„)* 



= 1 +e(.t), -1^x^+1 (65) 



in which x is a given constant, \x \ > 1, and the degree n of the dis- 
posable polynomial P(.t) is large. What is the approximate amplitude 
|e| of the equal error extrema of the minimax approximation? 

This is a special case of equations (55) and (56), which can be solved 
as a special case of equation (42) for which 

V(x) = ^ = (1 - x/x )K (66) 

To apply Section IV, define F(-) and D(-) by 

(1 - x/xtf = F(e iv ) + F(f iv ) = D(e iv )D(e- iv ), 

(67) 
F(e iv ) = 2C/'. 
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Log D(z) regular, \z\ ^ 1. 

We have already factored a linear function of £ in terms of exp (±i<p), 
in Section III. A similar factorization now gives 

W) = TTT^T . I 7 I < 1 (68) 

and then the coefficients of C a correspond to an expansion of 

a - 7 <na - ye-") 



1 +t" 



= ECe"* + EC/*'". (69) 

ff = ff = 



To determine e we only need the coefficients C a for a > n, which we have 
assumed to be large. 

The following expansion of [1 — y exp (up)], valid for |y| ^ 1, is well 
known 

(1 -ye iv )* = iLKJ'*; 

K = -1, K, = -t/2; (70) 

-(2.-3)! ^ 2 

4 ff_1 ((T - 2)!tr! 

When n is large and X « n, 

K n+ x +l = 2(w + X) - 1 
/C +x 2(n + X) + 2 

in which 



y ^ky (71) 



2n + 1 , 3 ,__. 

* 2n + 4 i 271 + 4- 1 " U ^ 

As a result, when n is large 

(1 - ye 4 *)* S £ /<V + ^''"T - < 73 > 

ff=n 1 — /eye v 

Now note that 

[ l - ye^Y \ K n+1 1 = A ,, (1 - fey')*g. +1 

L 1 + 7 2 J Ll - kye^J & ' "•" (1 + tV<1 - fc7e'")' 

(74) 

If equation (74) is used to evaluate C„ in equation (67), only the last 
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term contributes to C„ when a > n. Then 

t C +1+ xe' x * ~ (1-fcryg^ (75) 

fe (1 +7 2 )'(1 -W) 

which is a special case of equation (32) with A(-) a linear polynomial 
and B(-) a constant. The corresponding X(-) in equation (47) is a con- 
stant and cancels out. Then equations (68) and (75) applied to equation 
(50) give 

-2K. +1 



(1 -fcV)(l -ley 2 )*' (76) 

The constants Gi and G 2 contribute to the two highest degree terms in 
the polynomial P(-). They need not be computed unless the specific 
polynomial is needed as well as the amplitude |e| of the approximation 
errors. 

The linear A(-) and constant B(-) determined by equation (75) can 
be used in equation (40) to approximate the truncation error for the 
polynomial approximation defined by equation (66). The corresponding 
error in equation (65) can be found bj r dividing by (1 — x/x ) in . This 
gives 

_ max 1 gy | _ (1 - ky 2 ){\ + fc | y \ 

|e| ~ (1-|t|) { ° 

If k — 1, r = (1 + |t|) 2 < 4. Actually k < 1, but further analysis 
indicates that r will not be significantly >4 when y" is small. 

When |.t | — > °° , W(.r) — > 1, C n+1+ \/C n+l — > 0, and e T is dominated by a 
single Chebyshev poljTiomial (which has equal extrema). Consistent 
with this our 7 — > 0, then G, , G 2 — > in equation (76) and r — * 1 in 
equation (77). 

In equation (74), K„ +i can be determined by the formula for K, 
in equation (70). However, the following simpler approximate formula 
may be more useful: 

n+l 

Kn+1 - 2W*(n+l)(n+l/4)*' (?8) 

The error amounts to about 0.3 percent at n = 2 and about 0.04 percent 
at ?i = 6. The derivation is related to, but requires more than substitu- 
tion of Stirling's approximation for the factorials in equation (70). 
Fig. 3 illustrates computed errors t(x) and e r (a;) corresponding re- 
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Fig. 3 — Illustrating Example 1: (a) e(x) and (b) cr(x). 

spectively to our approximately equal ripple solution and truncation 
of the Chebyshev polynomial series. The constants ±e and re deter- 
mined by equations (74) and (75) are included for comparison. The 
computations started with 

x = 1.025, n = 10 

for which, as computed by equations (76), (77) and (78), 

7 = 0.8, k = 7/8, 

K n+1 ^ 0.0006881, e ^ 0.0040678, r ^ 3.74. 

5.2 Example 2 

Let 

P(:r.) = (1 - x/x ) i/2 + t(x), -1 ^ x £ +1 (79) 

in which the degree n of P(.r) is again large and x is again a given 
constant, \x \ > 1. 

Since the function F(x) is the same as in Example 1, equation (75) is 
again valid. Now, however, W(x) = 1 and hence Section 4.1 (on un- 
weighted polynomial approximations) is appropriate. Applying equation 
(75) to equation (37) gives 

(1 - ky 2 )^ 

(80) 



' ~ (1 + Y 2 )*(l " &¥) 



in which N(-) contributes only to the highest degree term in P(-). The 
error ratio r turns out to be 
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r = 



max e T 



= 1 + k I 7 I < 2. (81) 



5.3 Example 8* 

Let 

(1 - z/xo) 1/2 P(aO = 1 + e(aO, -lgJiSl (82) 

in which degree n of P(z) is again large and x is given, |z | > 1. 
In the equivalent weighted polynomial approximation 

F(x) = wfe = (1 " xA ° ri - (83) 

Proceeding as in Example 1, one now gets 

DV) = (1 ± -* 2) \ , 

(1 - 78")' 

V n />* - (1 + 7 2 )^ n+1 
2. c„ +1+xe _ (1 _ fc ^ )i(1 _ W>)J 

_ (27 1+ D!7 n+1 ^ 7" +1 

n+1 2 n+1 n! (» + 1)!~ [*-(n + 5/4)7' 
fc = 2n+4 =1 _ 1 



2n + 4 2n + 4 

Then equation (49) gives 

2if n+1 (l - Wf 
1 - fcV ' 



Equation (35) is now 

N 1 e~ iv 



(84) 



e ->AT( e --) = - Nie ~\ , (85) 

1 — 7e 

6 7(1 ~ *)(! + 7*)* 
iVl " 2 1 - fc 7 2 



— - E^-'". (86) 



1 — ye 
The first 2n + 1 terms in this series contribute to P(-) per equation 



♦The author has encountered this problem in connection with two different 
circuit theory studies, which will be described in other papers. 
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(36). 


The remainder 


can 


be summed, 


to get 






8(e 


1 - Ye-'* 



(87) 

Evaluating the error corresponding to simple truncation now gives 

_ max | e T | _ (1 - | y |)(1 - /c 2 7 2 ) 
r " |e| = (l-k\y\)(l-ky 2 ) (88) 

When n is large, k = 1 and r is so close to unity that the minimax 
refinement of simple truncation is not likely to be justified. However, 
our analysis has been useful in disclosing this fact, without the detailed 
computation of any minimax approximations. 

5.4 Example 4 

Previous work, which we shall discuss in Section VI, concerns the 
following problem: Let 

P n (x) = P B+ ,(aO + t(x), -1 ^ x £ +1 (89) 

in which P n+¥ (x) is a given polynomial of degree n + v and P n (rc) is a 
disposable polynomial of degree n. For what P„(x) does t{x) have the 
equal ripple form? 

Equations (32), (33) and (37) now simplify to 

ZX +1+x (e lx ") = IX + , + x(0 = B(e^), 

(90) 



l/CrJ _ xy*) 



X(e- iv ) 



+ «(e"), 



I Z(e <p ) + B(e ,> )I(e" > ) = e' i9 N(e~ ip ) 

in which #(•) is a polynomial of degree v — 1, with coefficients C n+1+ x 
and X(-) is a polynomial of degree i> — 1, to be found therefrom. The 
coefficients of B(-) can be found by expanding the left side of the last 
equation and equating to zero the coefficients of pgsitive powers of 
exp (i<p). The result can be expressed as the following specialization of 
equation (38): 

(c + | l)Q = (91) 

in which Q is again a column matrix whose elements are the v coefficients 
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of X(-) and / is the identity matrix of order p. The matrix C has the 
special form (assuming the elements q a of Q to be ordered per 
X(z) = £ q.z a ) 



C = 



'-'n+l '-'n + 2 

t^n + 2 ^n + 3 



c„ 



c B - 







c„, 












(92) 



When v = 1, X is a constant and the solution is elementary. When 
v = 2, X(z) is linear. Let z x be its zero. Then equations (91) and (92) 
require 



c 2 + 2C n+1 e - 4C 2 n+2 = 0, 



Z\ = 



2C n+i 



(93) 



The roots of the quadratic equation in e are real. When C n+i ^ 0, the 
larger \e\ > 2 | C n+2 | and | z x \ > 1, as required. Then 8 in equation 
(36) turns out to be a power series in exp ( — iip) which can be summed 
to get 

-i(2n + 2)»> 



*(0 = 



2B+2/7 



(94) 



1 - ye- ,ip 

7 = 1/2. • 

When C n+1 = 0, € = ±2C n+2 and ^ = ±1. But then the error due to 
simple truncation of the Chebyshev polynomial expansion of P„ +2 (.-c) 
is proportional to a single Chebyshev polynomial of degree n -f- 2, which 
has equal ripples with n + 3 extrema instead of n + 2. 

5.5 Example 5* 

As a last example consider the following nonalgebraic approximation: 



T(0) = £ A * sin o-^ = ^ + £(^) 

< 7T. 



(95) 



-ir < -e c ^ e ^ 

For what coefficients A, does the error e(0) have the equal ripple form 
and what is the amplitude e of the ripples? 



*This problem is of interest in, for example, the approximation of differentia- 
tion with a tapped delay line. A more detailed treatment, is planned for a future 
paper. 
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A sequence of transformations changes equation (95) into a weighted 
polynomial approximation. First, equation (95) is equivalent to 

n=l 

T{6) = sin 6 X) B p cos( P 0), (96) 

p = U 

2A P = B„_, - B p+1 . 

Second, relate to a new variable <p by 



sin 



- = q sin <p, 



< tt; (97) 



(98) 



q = sin -y < 1. 

Real ip maps into —d c ^6^d c . Also 

cos 6 = 1 — q 2 + <f cos 2<p, 

sin 6 = 2q{\ — q 2 sin 2 <p) h sin <p. 
In these terms, equation (96) becomes 

n-l 

p(0) = 2(/(l - q 2 sin 2 *>)* sin <* £ £„ cos (2pp) (99) 

0=0 

in which the set of coefficients B„ is linearly related to the set B„ of 
equation (96). In equivalent exponential terms 

t(6) = f (1 - q 2 sin 2 *) W) " P(«"")l ( 10 °) 

in which P{z) is a polynomial of degree 2n — 1 in odd powers of z only. 
Use equation (100) in equation (95) and solve for the factor in [ ]. 
The result can be expressed in terms of exponentials: 

P(«") - P(f i9 ) = F(e iv ) - F(e- iv ) 

+ D{e iv )D{e- iv )[E(e iv ) - B(« _< *)] (101) 

in which F(-), D(-), and E(-) are related to previous functions by 

- (1 - q 2 sin 2 <?)-> = £ 2C 2 ,_, sin (2<r - 1)*»; (102) 

W) = E^-/' 2 '" 1 "; 

i(l - g'sin 2 *) -1 = D(«")I>(«"'*); 
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r g g + tO 'P 



E(9) =±[E(e i n-E(e- i nV, 

It can be shown that the coefficients C 2ff -i obey a difference equation 
of order 2. The asymptotic behavior of the difference equation shows 
that 

bs a ->» , C 2 „ + , -» - (f^rf) 7<?2,-> . (103) 

As a result, for a sufficiently large n 



oo p i(2»+l)» 

if, 2 '- 16 = l+kye i2 * 

, ^ /2n + A 5 _ 4n + 1 



(104) 



2n 4- 3/ _ 4w + 3 



Proceeding almost exactly as in Example 3, using equation (104) and 
the -D(-) of equation (102), one can now obtain an approximation to the 
minim ax error e, to the error e T (<p) due to simply truncating the expan- 
sion of F(e" e ), and to the ratio r of max | e T | and | e |. As in Example 3, 
it turns out that the minimax approximation is only a little better than 
the approximation by truncation, at least when n is large. 

Figure 4 compares a computed e T (tp) with the approximate e and 
ratio r of max | e T | to | e |, using 

d c = 170° n = 15, 

for which 

7 = 0.8397 k = 0.96825, 

e^ 5.690° r^ 1.0840. 

VI. COMPARISON WITH OTHER WORK 

This section compares the present paper with previous publications 
in various related fields. It is not intended, however, to be an exhaustive 
survey of all related publications. 

The transformation from x to <p followed by distortion of the <p scale 
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TRANSFORMED f) SCALE 



Fig. 4 — Illustrating Example 5. 



to obtain an error function of the form e cos [(n + l)v» + /(<?)] has been 
used before. Pertinent references are papers by Clenshaw 2 and Stiefel, 3 " 4 
who call our j(tp) the "phase function". Of these, Clenshaw's paper is 
quite close to ours, and in fact our work might be regarded as a generali- 
zation of his. 

Clenshaw devotes much of his paper to the approximation of a 
polynomial of degree n + v with a polynomial of degree n, which is our 
Example 4. Clenshaw, (Ref. 2, pp. 30, 31) solves the problem for v = 2 
in a quite similar way, except that he expresses his Fourier series in 
terms of cosine functions instead of power series in e*'*. He exhibits an 
approximate solution which can be shown to be almost, but not quite 
equivalent to ours. We would have obtained an exact equivalent if we 
had restricted our polynomial P(-) to positive powers only, correspond- 
ing to a = to n in equation (25) instead of — n to n. In equation (36), 
this restricts the disposable G\& to X = 1 to n -f- 1 and increases the 
number of terms in 5(-) to X = n + 2 to *> . The result is a somewhat 
poorer, but frequently adequate, approximation to an equal ripple 
error. Clenshaw also notes how his approximation can be improved, but 
does not fill in the details. It can be shown that the improved approxima- 
tion would be an exact equivalent of ours. However, we have found that 
our formulations in terms of e ±lip , instead of cos <p, are simpler, and also 
more revealing concerning, for example, the nature of the approxima- 
tions. 

Clenshaw (Ref. 2, pp. 31-36) also considers v > 2, and obtains ap- 
proximate solutions for v = 3, 4 in terms of roots of cubic and quartic 
equations. However, he retains the use of cos <p, instead of e*'". As a 
result, he does not include a formulation for a general v in terms of an 
eigenvalue and eigen vector of a matrix, like our equation (91). 
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Clenshaw (Ref. 2, p. 29) notes that the equal ripple approximation to 
(x — x)' 1 has been solved exactly and cites Hornecker 6 and Rivlin 7 . 
Any solution to this problem is easily applied to more general problems 
in which the given function yields the same sort of remainder when the 
Chebyshev polynomial series is truncated. Examples are our Example 2 
and our general formulation for unweighted polynomial approximations 
with weight factor W = l,m= 1, /* = in the remainder function (32). 

Our procedures are more general in the following ways: First, remain- 
der functions can have the more general form (32). For practical pur- 
poses, degrees m and n should not be large. However, they need not be 
restricted to the special cases m = 0, m ^ 4 and m = 1, p =» 0. Second, 
minimax weighted errors can be obtained [by using suitable weight 
factors W(x) in equations (42), (44), and so on]. Third, unweighted mini- 
max approximations can be obtained with approximating functions of 
which the disposable polynomial is only a part (by solving an equivalent 
polynomial approximation with a weighted minimax error, as in our 
Examples 1, 3 and 5). Finally, relatively simple formulations have been 
obtained by using exponentials instead of cosine functions. 

Stiefel's papers 3 ' 4 have much less relevance to our work. They use 
the error formulation e cos [(n + l)<p + /(<?)] but obtain solutions by 
numerical iteration. Reference 3 also includes a general integral equa- 
tion, which determines the required coefficients implicitly but is not 
easily solved. 

Our use of rational fractions to approximate remainder functions, 
as in equation (32), and so on, is at least reminscent of the so-called 
e algorithm. The e algorithm also uses rational function approximations 
to remainders but for a different purpose — to increase the rate of con- 
vergence when functions are evaluated from their power series. It is 
quite different from the use of rational functions in the formulation of 
minimax polynomial approximations. References for the e algorithm are 
Shanks 8 and Wynn. 9 

Our procedures require evaluating certain of the coefficients in the 
Chebyshev polynomial expansion of a given F(x) (or in the equivalent 
Fourier series expansion in terms of <p). Various established numerical 
methods are available for this. 1 The best choice depends on the form in 
which F(x) is specified (for example in closed analytic form, as a power 
series in x, or numerically at a set of discrete points). When F(z) satisfies 
a differential equation with polynomial coefficients, the coefficients in 
the Chebyshev polynomial series are related by a difference equation of 
finite order and can be computed recursively. Our Example 5 is a special 
case. The general relation is described by Clenshaw 10 who also includes 
numerical tabulations of coefficients for some common functions. 
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The author's 1952 paper on network synthesis in terms of Chebyshev 
polynomials is only remotely related to the present work. 11 

VII. CONCLUSIONS 

Techniques like those described in Section IV and illustrated in Sec- 
tion V can be applied to many approximation problems in which the 
disposable part of the approximating function is a polynomial and 
approximately equal weighted or unweighted error extrema are desired. 
However, to be useful they must compete with other possible techniques, 
especially established numerical methods whereby equal-ripple ap- 
proximations are obtained by iterative improvement of a sequence of 
unequal-error approximations. This section notes some circumstances 
under which procedures like those described here may perhaps be pref- 
erable. 

First, the techniques described here are more likely to be competitive 
when the degree n of the disposable polynomial is large. When n is large 
iterative numerical methods are more likely to entail excessive amounts 
of computing. On the other hand, certain aspects of the more analytic 
techniques described here are likely to become easier as n becomes 
larger. These concern particularly the use of a simple rational fraction 
to approximate a remainder function, as in equation (32) . 

Second, the techniques described here are particularly suitable for 
exploring relationships between error amplitude | e |, the limits x a , x b 
of the approximation interval, the degree n of the disposable polynomial, 
and other parameters in the approximating function (such as x in 
examples 1, 2, and 3). In explorations of this sort the computation of 
the actual coefficients of the disposable polynomial P(x) can usually 
be omitted. When n is large this can mean omitting most of the computa- 
tions required for a complete determination of the approximating func- 
tion. Frequently, computations which end with | e | remain very simple 
even though n becomes arbitrarily large. 

Third, sometimes, as in our Examples 1, 2, 3 and 5, our techniques 
give quite simple estimates of the advantage of an equal-ripple approxi- 
mation over simple truncation of an infinite series of Chebyshev poly- 
nomials. Such a comparison may be useful, for example, in deciding 
what sort of approximation should be computed in detail. 

More generally, an attractive combination may be an initial explora- 
tion in terms of the techniques described here, followed by the detailed 
computation of one or more preferred cases by established iterative 
numerical methods. 

We have assumed here that the parameters disposable for purposes 
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of approximation are all the coefficients in a polynomial. Preliminary 
investigation indicates that similar methods may be feasible for dis- 
posable rational functions, or ratios of polynomials, provided the poly- 
nomials in the denominators are of quite modest degree. This will be 
the subject of a later paper. 
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